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ABSTRACT 


In  this  paper  Several  novel  aspects  of  radiation  cooling  are  discussed; 
particularly  those  associated  with  the  occurrence  of  surface  catalysed  atom 
recombination  at  high  Mach  numbers.  An  analogy  between  radiation  cooling 
and  chemical  surface  catalysis  is  explored  and  the  dominant  effects  resulting 
from  interactions  between  these  two  processes  are  illustrated  using  simple 
mathematical  models  which  serve  to  single  out  the  important  nondimensional 
parameters-  In  connection  with  flight  applications,  altitude « velocity  maps 
are  presented  which  provide  an  overall  picture  of  the  regimes  in  which  chemicai 
nonequilibrium  effects  should  be  noticeable.  Several  unexplored  areas  of 
potential  interest  in  the  design  of  hypersonic  lifting  vehicles  are  outlined  and 
related  to  the  available  literature.  Interestingly  enough,  not  all  of  these  areas 
deal  udth  thermochemical  effects.  Many  "classical"  problems  remain  unsolved, 
as  illustrated  by  a  brief  discussion  of  the  radiation  cooled  flat  plate* 
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1.  INTRODUCTION 

The  surface  temperature  at  each  point  along  an  object  traveling  at  high  speed 
throu^  a  planetarv  atmosphere  will  seek  a  value  at  which  the  rate  of  heat  input  by 
convective  processes  exactly  balances  that  lost  by  phase  change,  radiation  and 
internal  heat  conduction,  oepending  largely  upon  velocity  and  free  stream  density 
level,  thezmai  radiation  from  the  surface  may  be  adequate  to  prevent  appreciable  phase 
change  and,  can,  in  fact  become  the  dominant  form  of  heat  rmnovaL  It  is  therefore 
of  interest  to  inquire:  (i  )  what  local  surface  temperatures  would  result?,  and 
( ii  )  how  would  the  surface  temperatures  vary  i^th  environmental  conditions?  To  cite 

a  particular  example  of  current  interest  (l  -  3),  it  is  known  that  local  deceleration 

of  the  gaseous  medium  relative  to  a  high  speed  entry  object  gives  rise  to  temperature 

levels  at  v#iich  even  the  most  stable  molecules  are  fragmented,  (eg*  nitrogen  in  the 
aOnospheres  of  the  earth  and  mars),  this  may  occur  either  behind  the  bow  shock  wave 
that  Stands  off  blunt  nosed  bodies  or  in  the  Shear  layer  near  lateral  surfaces  as  a 
result  of  locally  intense  viscous  heating.  When  the  gaseous  medium  experiences  this 
kind  of  thermochemical  change  steady  state  surface  temperatures  may  be  expected  to 
depart  considerably  from  those  corresponding  to  the  aerodynamically  equivalent  motion 
through  a  chemically  "inert"  aOnosf^ere  (^).^  Because  of  the  dominant  role  which 
surface  temperature  plays  in  the  determination  of  physical  and  chemical  integrity* 


gas  cap  radiation  ^  the  surface  will  not  be  included  in  the  present  discussion. 

As  will  be  seen,  this  radiation  is  usually  negligible  in  the  flight  regimes  in 
which  cooling  solely  by  radiation  becomes  feasible, 

^an  "inert"  atmosphere  can  be  visualized  as  one  in  which  the  kinetics  of  the  thenno* 
chmnical  changes  (e,g.  dissociation,  ionization)  are  not  sufficiently  fast  for  object* 
in  the  size  range  of  interest. 
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particular  ern^asis  must  be  placed  on  its  accurate  prediction  using  available  physico-* 
chemical  infosnation,  experimental  heat  transfer  data  (usually  obtained  in  the  absence 
of  complete  simulation)  andj  of  coursei  theory^  However  our  purpose  here  is  not  so 
much  to  present  accurate  GOffiputational  methods  as  it  is  to  Illustrate,  by  means  of 
relatively  simple  examples,  what  dominant  effects  are  likely  to  be  encountered  and 
in  What  flight  domains.  This  admittedly  limited  goal  is  nowhere  more  evident  than  in 
our  treatment  of  thermo  chemical  effects,  where  it  v41l  be  appreciated  that  the  problem 
of  making  accurate  predictions  is  computationally,  if  not  conceptually,  complex. 

While  the  radiation  cooling  of  lifting  re-entry  vehicles,  mierometeorites,  and  jet 
aircraft  is  implicit  in  much  of  the  discussion  which  follows  it  should  be  mentioned 
that  similar  considerations  would  apply  to  the  radiation  cooling  of  nuclear  rocket 
nozzles  in  Outer  space  (5)  or  in  more  mundane  applications,  such  as  the  determination 
of  radiation  errors  in  the  pyrometry  of  flame  gases  (6),  or  the  laboratory  use  of 
differential  catalytic  probes  in  the  deteimination  of  local  free-radical  concentrations 
in  npneQui librium  flow  systems  (2,  §)• 

Happily,  many  important  aspects  of  radiation  cooling  are  now  well  known, 

(9  -  1^)  and  graphs  have  become  available  to  facilitate  the  deteimination  of  steady 
state  surface  temperatures  at  various  points  along  aerodynamic  shapes  of  simple  geo^ 
metric  foim  (U  -  13).  However  the  design  problem  is  by  no  means  solved  since 
available  results  are  far  from  universally  applicable.  Indeed,  technological  advance 
continues  to  introduce  new  parameters  at  a  rate  which  poses  an  intellectual  ehillenge 
which  must  be  met  if  this  advance  is  to  be  sustained  at  reasonable  cost. 

For  definiteness,  consider  the  surface  temperature  predictions  shown  in  Fig,  1, 
constructed  from  the  recent  work  of  Naysmith  and  Woodley  Here  are  epllected 

the  predicted  steady  state  temperatures  of  slender  surfaces  whose  emissivity  is 
taken  as  0,3  (being  representative  of  polished  metals).  Results  are  applicable  to 
Streamwise  location  x  =  I  ft.  and  cover  the  flight  Mach  number  range  5  to  10  at  twm 
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altitudes^  Three  important  features  are  immediately  evident.  First*  surface  tern*' 
peratures  are  considerably  below  the  corresponding  gas  dynamic  recovery  temperatures, 
particularly  for  high  altitude  flight.  Thus,  for  Mach  10  flight  at  200  Kft,  the 
local  surface  temperature  is  estimated  to  be  only  S40®k  whereas  the  gas-^dynamic 
recovery  temperature  in  the  absence  of  dissociation  would  be  about  6060®K  at  this 
same  flight  condition.  If  One  were  limited  by  a  sudden  loss  of  strength  or,  say,  by 
local  melting,  neglect  of  radiation  would  therefore  provide  a  very  pessimistic  esti¬ 
mate  of  the  maximum  attainable  Mach  number  based  on  aerodynamic  heating  considerations 
alone,  second,  in  the  absence  of  internal  heat  eonduction  considerable  surface  tem¬ 
perature  nonuniformity  would  exist,  as  evidenGed  by  the  fact  that  the  local  temperature 
is  always  much  lower  than  the  (mean)  temperature  which  would  have  existed  had  the 
total  heat  imput  from  the  leading  edge  to  the  point  in  question  been  distributed 
uniformly  over  that  distance.  For  the  case  cited  above  this  temperature  difference 
amounts  to  Some  160°K,  The  very  existence  of  large  streamwise  temperature  variations, 
especially  for  laminar  boundary  layers,  raises  the  question  of  the  validity  of  using 
isothermal  surface  heat  transfer  GoefflGients  (^).  This  question  will  be  discussed 
in  Section  2*2.  Finally,  surface  temperatures  are  very  sensitive  to  the  fluid  dynamic 
state  of  the  boundary  layer  and  are  much  larger  if  the  boundary  layers  are  turbulent, 
reflecting  the  ertianced  conditions  of  heat  input.  At  M^^  =  10,  h  =  50  Kft,,  for 
example,  a  turbulent  boundary  layer  would  eause  ipcal  Surface  temperatures  to  be  some 
490°K  larger  than  in  the  case  of  laminar  boundary  layers,  bringing  the  steady  state 
surface  temperatures  to  above  2000°K. 

While  the  error  in  these  results  was  estimated  by  the  authors  to  be  only  about 
10  percent,  at  intermediate  altitudes  and  high  Mach  numbers  the  kinetics  of  oxygen 
dissociation  and  recombination  could  begin  to  influence  results  of  this  type  consi¬ 
derably,  as  would  be  suggested,  for  example,  by  the  difference  between  the  real  and 
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ideal  stagnation  tempefatures  computed  at  the  free  stream  pressure.  In  fact,  it  will 
he  shown  that  the  steady  state  surface  temperature  can  be  either  higher  or  lower  than 
that  estimated  on  the  basis  of  ideJl  gas  calculations  (l^).  Moreover,  owing  to  the 
possibility  of  exothermic  atom  recombination  within  the  boundary  layer  or  at  the 
radiating  surface  itself  the  surface  temperatuTe  at  a  nose  cap  or  leading  edge  can 
certainly  be  much  greater  than  the  real  gas  stagnation  temperature.  ThuSj  much  of 
the  data  contained  in  available  design  charts  should  be  interpreted  in  the  light  of 
■Uie  chemical  kinetic  considerations  (ef.  Part  3),  Finally^  in  realistic  situations 
chemical  kinetics  can  introduce  entirely  new  phenomena  [e^gs  surface  temperature 
hysteresis  (j^)]  as  well  as  modify  the  "expected"  dependence  of  surface  temperature 
on,  say,  physical  scale  New  design  parameters  (e.g.  recombination  coefficients) 

are  thereby  introduced  into  the  materials  selection  problem,  and,  as  such,  this  class 
of  phenomena  merits  further  attention. 

2.  RADIATION  COOLING  IN  THi  AiSINCE  OF  CHEMICAL  REACTION 


In  a  remarkable  paper  on  the  theory  of  laminar  boundary  layers,  Lighthill  (g) 
pointed  out  that  since  the  wall  temperature  distribution  T^(x)  of  a  radiation  cooled 
body  is  itself  non*uniform,  a  theory  strictly  applicable  under  this  nonisothermal 
condition  is  required  to  compute  it.  Using  a  method  which  is  in  fact  asymptotically 
exact  in  the  limit  of  infinite  Prandtl  number^  he  was  able  to  show  that  the  true  steady 
state  temperature  distribution  is  obtainable  from  the  solution  of  a  nonlinear, 

singular  Volterra  integral  equation.  For  a  flat  plate  of  zero  thermal  conductivity 
and  constant  emissivity,  slowly  convergant  expansions  for  the  steady  state  temperature 


for  slender  objects;  for  blunted  shapes  this  difference  is  reduced  owing  to  the 
higher  pressure  levels  encountered, 

^While  Lighthill  discussed  applications  to  laminar  flow  only  in  the  limit  Pr^  ^  ^ 
his  basic  solution  is  applicable  to  turbulent  flow  as  well,  since  the  thermal  boundary 
is  then  fully  immersed  in  the  laminar  sublayer  (16). 
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distribution  were  obtained  (for  small  x  and  large  x)  containing  fraGtional  powers  of 
the  streamwise  coordinate  x.  Lighthill  did  not  investigate  the  overall  radiative 
transfer  from  such  plates,  or  the  validy  of  aipproxifflate  methods  Miiieh  make  looai  use 
of  ''isothermal"  heat  transfer  coefficients,  Nuj^gg(x)i  Yet  it  is  clear  from  Lighthilr'S 
Series  solution  that  near  the  leading  edge  of  the  plate  Nu/Nu^^^  approaches  the  finite 
value  1436^,  with  Nu/nuj^^^  approaching  unity  only  far  downstream  (as  x^  00  )i  in 
195®  N.  Curie  (^)  applied  an  approximate  method  for  calculating  the  surface  temper* 
ature  distribution  of  a  radiation  cooled  flat  plate »  While  his  approximation  reduces 
the  problem  to  an  algebraic  (4th  degree)  equation,  it  is  readily  shown  to  be  far  less 
accurate  than  straightforward  application  of  NUj^^^(x)  to  the  same  problem;  a  technique 
which  also  leads  directly  to  an  algebraic  equation  of  identical  focmi  In  the  light 
of  this  background  the  remainder  of  this  section  will  therefore  be  devoted  to  (i)  a 
brief  examination  of  the  accuracy  of  the  NUj^^(x)  approximation  for  both  laminar  and 
turbulent  boundary  layer  (ii)  the  development  of  an  improved  approximate  method  readily 
applicable  to  laminar  or  turbulent  boundary  layer  flow  over  surfaces  with  temperature 
dependent  emissivity,  and  (iii)  a  comparison  between  the  overall  radiative  transfer 
from  plates  of  zero  and  infinite  theimal  conductivity  (lO) .  It  will  be  shown  that 
the  radiation  cooling  problem  is,  mathematically,  strictly  analogous  to  a  problem  in 

the  diffusional  theory  of  surface  catalyzed  chemical  reactions  (w,  18,  19).  This 

analogy  will,  in  fact,  provide  a  natural  link  with  the  remainder  of  this  paper,  devoted 

in  part  to  the  effect  of  exothertnie  surface  reactions  in  radiation  cooling  situations. 

Z,\  A  Further  Analogy  in  Transport  Theory  (lo) 

Following  Lighthill  (f),  let  us  begin  our  discussion  of  radiation  cooling  at 
hi^  Mach  numbers  with  a  treatment  of  the  incompressible,  M„^  0,  planar  problem' 
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In  that  case  the  radiation«cooiing  chemical  surface  catalysis  analogy  is  readily 

demonstrated  as  follows.  OGHsider  the  flow  configurations  depicted  in  panels  1  and 

2  of  Fig.  2.  Panel  1  shows  the  incompressibie  flow  of  a  hoti  nonabsorbing  fluid 

* 

over  a  convex  cylinder  of  zero  theitnal  conductivity.  In  the  steady  state  each 
element  of  surface  will  radiate  heat  away  at  exactly  the  rate  reaGhing  it  by  convec¬ 
tion  from  the  main  streami  i.e. 


(U 


Here  x  is  the  thermal  conductivity  o^f  the_ fluid ,  t  the  total  hemispheric 
emissivity  of  the  surfacOj  ^  is  the  stefan-Boitzmann  raidiation  eonstantj  and  T 
represents  the  local  fluid  temperature  (assumed  equal  to  the  solid  surface  temperature 
at  y  «o).  In  general,  &  is  itself  tmnperature  dependent  so  that  the  product 

W 

^ taken  to  vary  as,  say,  the  nth  power  of  local  surface  temperature, 
where 

It-  M  ^  ^  cl  S ^  (J 


For  simplicity,  it  will  be  assumed  herafter  that  n  may  be  taken  to  be  constant  over 
the  surface  temperature  range  of  interest,  although  different  from  4.  Some  typical 
values  of  n  are  collected  in  Table  1,  where  it  is  seen  that  values  between  2  and  6 
are  conceivably  of  practical  interest,  with  most  ceramic  materials  exhibiting 
n  <  4  and  metals  exhibiting  n  >  4, 


^  .... 
m  practice,  this  condition  may  be  replaced  by  the  condition  of  sufficiently  thin 

walls  or  the  use  of  segmented  walls  to  prevent  heat  conduction  in  the  streamwise 

direction. 
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At  sufficiently  large  Reynolds  numlDerS  the  velocity  field  in  the  neighboriiood 
of  the  Cylinder  will  satisfy  the  constant  property  boundary  layer  equation 


U 


'2>il 


du. 


Hx 


71C. 


(3) 


where  i)  and  6^^  are^  respectivelyj  the  kinematic  viscosity  and  eddy  viscosity  of 

the  fluid  and  u^(x)  is  the  main  stream  velocity  distribution  around  the  cylinder. 

Similarlyi  at  sufficiently  large  Peelet  numberSj  the  tempefature  field  T(xiy)  vd,thin 
the  fluid  will  satisfy  the  constant  property  boundary  layer  equation 


?T  ^  r 

-h  y  ^  ^  ^  , 

'dl  0(j. 


9^/  I  rr  ‘‘‘  li)  J  'i  it 

x\  ,\  d 

represent  the  laminar  and  turbulent  Prandtl  numbers  for  heat 
conduction*  Now  consider  panel  2  of  Fig.  2*  Depicted  here  is  the  steady  flow  of  an 
incompressible  carrier  fluid  over  an  impeimeable  catalytic  cylinder.  Contained  in 
the  carrier  fluid  is  a  reactant  (present  in  dilute  amounts)  v^ich  diffuses  to  the 
solid  surface  and  reacts  there  at  a  rate  which  is  some  function  of  the  local  reactant 
coneentrltiOn  (mass  fraction)  (j^  itself.  Very  frequently  this  phenomenological 
reaction  rate  law  is  well  represented  by  a  simple  power  law,^  so  that  the  reactant 
concentration  field  at  the  fluid/solid  interface  must  satisfy  the  boundary  condition 


(4) 


where  Prj^  and  Pr^^ 


(t) 


(5) 


* 

The  validity  of  the  analogy  described  here  is  not  contingent  upon  the  validity  of 
boundary  layer  theory}  however*  the  analogy  will  be  Illustrated  later  for  the  case 
of  boundary  layer  flow  over  a  flat  plate. 

*an  important  example  for  which  this  is  not  true  is  given  in  Section  3.3  (23) 
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called  the  ihtr^insic^  catalytic  activity  of  the  surface^  h  is  the  true 
rgaGtion  otdef  and  D  is  the  Fick  diffusivity  for  reaGtant  diffusion  through  the 
carrier  fluid.  For  a  Single  material  at  uniform  surface  temperature  the  two  para^- 
meters  and  n  then  Gompletely  define  the  chemical  kinetic  characteristics  of  the 
Surface. 


The  velocity  field  in  the  neighborhood  of  the  catalyst  will  be  identical  to  that 
in  panel  1  and,  moreover,  in  the  absence  of  ehemical  reaction  in  the  fluid  i^ase  the 
conGentration  field  will  satisfy  the  constant  property  boundary  layer  eguation 


9  a' 
u.  t: 


+  1/ 


f 


^C\ 


f 


I 


f  ^ 

-r  -h 


p 


\i^  1 


(6) 


Far  from  the  catalyst  (y  oo  )  we  have  c^(x,y)  U  s  const.  The  mathematical 

e 

similarity  between  the  two  problems  is  then  evident  Since  T(xjy)  and  o6(x,y)  Satisfy 
differential  equations  and  boundary  conditions  of  identical  form.  In  particulari 


T^(x)  H  T(x,o)  and  OC^M  i  0^(x,o)  are  not  known  a  pri^i  but  in  fact  constitute  an 

important  part  of  the  solution.  Thus,  having  Obtained  T  (x),  one  Can  calculate  the 

w 

local  and  overall  rates  of  heat  transfer  (per  unit  depth)  from  the  relations 


f'fx)  «  (X)  /  ^ 


(7) 


Q  f  (T' Jz 

'  o 

(for  one  side  of  a  symmetrical  cylinder).  Similarly,  having  obtained  dne 

W 

can  obtain  the  local  and  overall  rates  of  reaction  (per  unit  depth)  viz. 


9 


TP-60 


^"(x)  = 

and 

o 

(for  one  side  of  a  syntmetricai  Gylinder).  The  exponent  n  defined  by  iq.  (2)  plays 
the  same  role  in  the  radiation  cooling  problem  as  the  true  reaGtion  order  n  in  the 
diffusion  surface  catalysis  problem^  Interestingly  enough,  v4iile  the  range  2  to 
6  was  seen  to  be  of  interest  for  radiating  surfaces  (cfi  Table  1)  in  heterogeneous 
chemical  kinetics  one  usually  encounters  values  of  the  true  reaction  order  n  between 
0  and  2i  Thus,  solutions  covering  the  range  o  <  n  <  6  would  be  of  general  interest 
in  transport  theory.  When  expressed  in  terms  of  nondimensional  transfer  rates  and 
coordinates,  the  solution  for  n  =  2  would  then  be  directly  applicable  to  either 
problem.  The  nature  of  the  important  nondimensional  variables  will  be  evident  in 
discussing  a  particular  case,  viz,  the  radiation  cooled  flat  plate. 

2.2  Approximate  Solution  for  the  Radiation  Cooled  Flat  Plate?  Incompressible  Case  (10) 
An  approximate  method  which  has  already  proven  convenient  for  solving  the 
diffusion-surface  reaction  problem  will  be  outlined  here  and  then  used  as  the  basis 
for  a  brief  discussion  of  accuracy. 

Ambrok  (24)  has  provided  a  simple  expressiPf'  fpr  estimating  the  effect  of 
surface  temperature  variations  T  (x)  [or,  more  generally,  ^(x)  i  T_  -  T„3  on  the 
local  heat  transfer  coefficient.  For  a  flat  plate  (cf.  Fig,  3)  his  result  can  be 
cast  in  the  Stanton  number  fopn 
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J(xU)l 


where  the  eXiponefit  ^  takes  eh  the  value  2  f&t  laminar  boundary  layer  flOWj  and  5 

for  turbulent  boundary  layer  flow  (tripped  at  leading  edge)»  St^g^(L)  represents 

the  appropriate  value  of  the  Stanton  number  at  the  trailing  edge  (x  =  l)  in  the 

event  Si  =  const »  The  success  of  this  method  can  be  conveniently  tested  against 

exact  (similar)  solutions  for  the  case  ^  ^  ^  where  1  is  constanti  Results 

over  the  1* range  of  interest  heie(o  <  l  <  j)  are  shown  in  Fig»  4.  Svidently,  the 

method  is  quite  satisfactory  for  both  torbuient  and  laminar  boundary  layer  flow  pro- 

.  * 

vided  the  Prandtl  number  is  not  too  small*  This  being  the  easej  consider  the 
application  of  Eq,  (ll)  to  the  radiation  cooling  problem.  Using  the  definition  of 
St(x)  the  boundary  condition  (1)  takes  the  fo«n 


where  St(x)  is  given  by  Iq.  (ll)  with  cd-  2  or  5.  This  nonlinear  integral  equation 
can»  however»  be  converted  to  a  separable  ordinary  differential  equation  and>  thus* 
readily  solved.  At  this  point  it  is  convenient  to  introduce  the  following  nondi* 
mensional  variables 


Sn  the  laminar  case  it  is  well  known  that  the  Pr  ^  op  asymptote  provides  an  accurate 
estimate  of  translBr  coefficients  even  for  Prandtl  numbers  of  order  unity  (|^. 
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^  '  CT"  /T 


4 


.  rar/z  ; 


(U) 


Note  that  X  is  a  stretched  streamwise  coordinate  proportional  to  the  local  boundary 
layer  thickness  terms  of  these  variables  (12)  becomes 


^  M  Z 


OJ^i 


OJ  =‘/ 


LU 


which,  upon  differentiation^  leads  to  a  Separable  differential  equationi.  This 
equation  can  then  be  integrated  tertn-by-term  as  follows 

to  jf 


Co 


CO 


<JU  -/ 


^Cu  yt 


J  CO -f  ^  ^  j  -  CO 

6 

giving  the  steady  State  temperature  distribution  in  the  inverse  form  z  #2(^5  n,  Ct? ) 
for  any  value  of  n  and  for  2,5,  lesults  are  shown  in  Fig.  5  for  the  cases 

n  =  2, 3,4, 5,6.  Solid  curves  pertain  to  laminar  boundary  layer  flow,  whereas  dashed 


contours  correspond  to  turbulent  boundary  layer  flow.  It  should  be  remarked  that  if 


one  had  neglected  the  effect  of  T  (x)  on  the  local  heat  transfer  coefficient  Nu(x) 

W 

then  laminar  and  turbulent  results  would  be  identical  in  theee^eooLdi nates .  While 
the  difference  does  not  appear  too  impressive  in  Fig.  5,  it  is  amplified  in  calcu-* 
lating  the  local  rate  of  radiative  transfer  (normalized  by  the  value  at  the  leading 
edge),  i.e.  q  1  @  this  is  shown  in  Fig,  6  for  the  special  case  n  ^  4  (constant 
emlssivity).  Laminar  boundary  layer  results  for  the  case  n  -  4  compare  very  favorably 
with  the  exact  solution  (for  Pr=^  00  )  given  by  Lighthlll  (2),  with  departures  not 
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exceeding  0.17  pereent  over  the  limited  range  of  z  Govefed  by  Lighthill's  p6wer*sefies 
expansions.  As  suggested  fay  Fig.  4«  still  better  agreement  would  be  expeGted  in  the 
turbulent  Gase^  sinoe  heat  transfer  Goefficients  are  then  less  sensitive  to  surfaGe 
temperature  variations. 


of  partiGular  interest  is  the  overall  rate  of  radiative  transfer  for  the  entire 
surfaGe  [of.  Iq.  (l^)].  This  is  conveniently  desGribed  in  terms  of  a  nondimensional 
GoeffiGientt  t)>  defined  as  follows 

<3- 

This  quantity  has  the  eharaeter  of  a  radiative  "fin  effiGienGy**  (26).  Mien  expressed 
in  terms  of  the  dimensionless  variables  ^  *  z»^  one  finds 

Ji 

But  an  expression  for  qiz'^'^^dz  in  terms  of  Q  alone  is  readily  obtained  from  Eq.  (17) 
So  that  teim-by>tetm  integration  is  possible.  One  thus  obtains  the  remarkably  simple 


result 


(20) 


wbere*  in  this  ease»  B  and  z  are  evaluated  at  the  trailing  edge  {z  =  J\j  ).  Results 
for  the  integrated  radiation  effectiveness  factor  are  collected  in  Fig.  S  where  the 


absieca  is  not 


but  rather 


(21) 
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With  this  choice  of  afesiccai  it  is  seen  that  laniinaf  and  turbulent  results  are  almost 
GoinGident  for  each  value  of  the  parameter  n.  The  parameter  is  merely  a  dimen¬ 
sionless  radiation  parameter  identical  to  ^  [cfi  iq.  (14)3  ^bt  based  on  the  global 
heat  transfer  coefficient  Nu.  ^  rather  than  Nu.  (L). 

ISO  ISO 

we  are  now  in  a  position  to  make  an  interesting  comparison  between  two  extreme 
caseSi  Imagine  two  flat  plates  differing  only  in  their  thermal  conductivities  and 

placed  in  identical  fluid  dynamic  environments  *  in  the  steady  state*  which  plate  will 

lose  more  energy  by  radiation,  and  thereby  extract  energy  at  a  greater  rate  from  the 
fluid?  Ividently,  this  comparison  is  most  readily  carried  out  for  the  extremes  of 
zero  and  infinite  thermal  conductivity^  In  the  latter  case,  by  an  overall  energy 
balance,  the  perfect  conductor  must  assume  a  uniform  temperature  given  by  the  solution 
of  the  algebraic  equation 


The  relevant  radiation  effectiveness  factor  q  would  then  be  computed  from  ©^.  When 
this  comparison  is  made  one  obtains  a  very  intriguing  result;  viz,  while  the  insulating 
plate  always  radiates  away  more  energy  than  the  conducting  plate,  the  difference  is 

extremely  small  for  all  values  of  the  radiation  parameter  A.  For  constant  emis- 
Sivity  flat  plates  the  difference  never  exceeds  2^  percent  in  the  laminar  case 
(attained  at  ^  ^  1.2)  or  0.27  percent  in  the  turbulent  ease  (attained  at^K  1,5), 
Moreover,  this  difference  will  decrease  as  n  increases  so  that,  for  practical  purposes, 

it  may  be  said  that  metals  of  either  zero  or  infinite  thermal  conductivity  would 
radiate  away  equal  amounts  of  energy  per  unit  time. 

Of  course*  in  practice*  a  very  important  difference  in  these  two  cases  is  the 
existence  of  very  high  surface  temperatures,  and  gradients,  in  the  leading  edge  region 
of  insulating  surfaces.  This  may,  for  instance*  cause  local  melting  of  a  low  conduc¬ 
tivity  plate  under  conditions  for  which  a  plate  of  much  larger  thennal  conductivity 

would  remain  intact.  However*  two  comments  should  be  made  about  the  leading 
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edge  fegidn  in  view  of  the  aipproximatione  inherent  in  the  theory.  Firsts  as  pointed 
Out  by  Lighthill  (g),  the  infinite  gradients  appearing  atrx  =0  are  a  GOnsequenGe 
of  the  predicted  behavior  of  the  boundary  layer  thickness  ^  ^  This 

behavior  is  incorrectly  given  by  boundary  layer  theory  so  that  the  surface  teniperatures 

&(z}n)  indiGated  in  Fig.  5  for  z  ^  d  will  in  fact  be  rounded  off  at  a  somewhat  lower 

a 

value.  Liglithlll  estimated  this  effect  by  noting  that  the  boundary  layer  at  x  =  o 
can  be  regarded  as  being  of  finite  thickness^  or  order  9/o  for  Pr  =  0(1).  Balane* 
ing  the  heat  input  X[t  -  T  (+o)]iu  /-p  against  the  radiation  loss  (+d)]* 

gives  an  equation  for  d(d)  of  the  form  (22)  but  with  ^  replaced  by 


(23) 


Consequentlyt  when  this  parameter  is  smalU  we  have 

e(o)  i  - 

Secondj  owing  to  the  large  gradients  occurring  near  the  leading  edge,  longitudinal 
heat  conduction  would  be  expected  to  play. Some  role  in  reducing  the  actual  temper* 
ature  in  this  region.  But>  since  the  nondimensional  parameter 


% 


L 


ISO 


is  frequently  very  small  (vdiere  is  the  thermal  conductivity  of  the  solid  and  t  is 


its  thickness  in  the  leading  edge  region),  this  suggests  that  the  effects  of  finite 
material  conductivity  will  be  confined  to  a  leading  edge  "boundary  layer"  (within 


the  material)  of  axial  extent 


*Of  course*  for  blunted  shapes  finite  so  that  @(o)  can  be  considerably 

less  than  unity  (cf.  Section  3,3). 
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L 

Thus*  conduction  alone  would  be  exj^cted  to  cause 


to  differ  from  unity  by  an 


amount  proportional  to  the  small  parameter  given  in  ipi  (25) ^  Returning  to  the 
radiation  cooling/catalytiG  surface  reaction  analogy*  the  first  of  these  effects 
(i.e^  breakdown  of  boundary  layer  theory)  will  also  be  applicable  for  catalysts  with 
sharp  leading  edges;  however*  the  correction  for  axial  diffusion  along  the  solid  is 
negligible.  With  these  restrictions  in  mind  we  may  then  conclude  this  section  wdth 
a  Synopsis  of  the  analogy  and  a  brief  statement  of  the  accuracy  of  computing  T^(^) 
using  local  values  of  the  isothermal  heat  transfer  coefficient  Nu.„  (x). 


Table  2  sutnnarizes  the  radiation  coo ling/cata lytic  surface  reaction  analogy  for 
laminar  or  fully  turbulent  boundary  layers  developing  along  a  flat  plate.  It  is 
seen  that  the  analog  of  ^  in  the  catalysis  problem  is  the  nondimensional  catalytic 

parameter  (22>  ^) 

which  compares  the  characteristic  rate  of  reaction  (numerator)  to  the  characteristic 
rate  of  reactant  diffusion  (denominator).  In  the  diffusion* reaction  case  the  Nusselt 
number  i§  numerically  egual  to  the  heat  transfer  coefficient  but  can  be 

obtained  from  it  by  making  the  parameter  replacement?  Pr^  =►  Pr^.  Here  Prp(i  ^/p) 
is  the  diffusional  Prandtl  number.  The  nondimensional  transfer  coefficients  a  and 
a  introduced  in  Table  2  may  be  used  in  place  of  t|*  t)  in  some  applications.  For 
instance*  I  compares  tiie  actual  rate  of  radiative  transfer  with  the  total  convected 


*In  the  Western  chemical  engineering  literature  usually  called  the 

Sherwood  number  and  Pr^  is  called  the  Schmidt  number. 
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fate  at  which  heat  would  have  to  be  removed  to  maintaih  the  plate  at  a  tempefatufe 
Mhich  is  neglifible  by  Gompafisoh  with  The  Goefficients  fj  and  5  are  cleatly 
felatedj  in  fact 


It  may  be  verified  that  eaGh  of  the  curves  shown  in  Figs  7  has  the  pfOpefty 
for  large  Showing  thatj  in  this  extreme,  5  =*■  1*  Againi  by  analogy  with  the  diffu» 
sion-surface  reaction  case,  it  is  conceivable  that  the  theoretical  relations  discussed 
in  this  section  could  be  used  to  experimentally  infer  total  hemispheriC“emissivitieSi 
Such  an  application  would  be  akin  to  the  measurement  technique  described  in  iefi  (^) 
where  the  steady  state  temperature  distribution  in  a  radiating  fin  is  used  to  obtain 
its  surface  emissivityi 

While  it  can  be  shown  that  the  local  application  of  NUj^^(x)  will  always 
underestimate  the  actual  local  surface  temperature,  the  magnitude  of  this  effect  is 
relatively  small,  especially  for  turbulent  boundary  layers  with  n  ^  4  metals). 
However,  for  laminar  boundary  layers  developing  along  radiating  solids  with  g  «  2 
(ceramics)  the  errors  can  reach  six  percent  for  Pr  ^  0(1)  and  exceed  this  value 
for  small  Prandtl  numbers.  It  can  also  be  shown  (30)  that  the  errors  incurred  for 
flat  plates  cannot  be  exceeded  in  the  presence  of  non«zero  pressure  gradients  (e.g. 
wedge  flows).  Therefore,  inmost  practical  applications,  local  application  of 
isothermal  heat  transfer  coefficients  should  result  in  adequate  surface  temperature 
predictions, 

2.3  Extension  to  High  Mach  l*imbers  (9) 

For  the  case  of  laminar  flow  over  a  flat  plate  at  low  Mach  numbers,  viscous 
dissipation  of  kinetic  energy  begins  to  modify  the  steady  state  temperature  profile 

and,  hence,  the  total  rate  of  radiation  loss. 
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before  significantliy  altering  the  fluid  properties.  POlhausen  (31)  showed  that,  In 
ths  abse'Ocs^  of, heat  Aransi.er,  the  surface  temperature  would  then  achieve  a  uniform 
value,  T_,  which  is,  in  general,  higher  than  Tg  and  different  from  the  mainstream 
total  temperature 


(m 


By  numerically  solving  a  generalized  form  of  Eq.  (4),  viz . 


9r 


the  recovery  temperature  was  found  to  be  given  by  an  expression  of  the  form 

2. 


T  ^  r  +  r 

f  ^  -  I  f  f  ^ 


) 


4-  ii 


(31) 


vAtere  the  recovery  factor  r  ^(P\)  represented  by  f®*"  values  of  the 


Prandtl  number  not  very  different  from  unity.  Owing  to  the  linearity  of  Eq*  (30) » 
solutions  of  the  homogeneous  energy  equation  (4)  can  be  added  to  this  special  solution 


of  the  inhomogeneous  equation  (30)  and  still  satisfy  Eq,  (30).  Thus,  heat  transfer 

in  the  presence  of  viscous  dissipation  can  be  calculated  using  =  0  heat  transfer 

coefficients  by  merely  replacing  the  mainstream  fluid  temperature,  Tg,  with  the 

recovery  temperature  Tj.[cf.  Eq.  (31)3*  At  still  larger  Mach  numbers  the  theimal 

properties  of  the  fluid  can  no  longer  be  considered  constant,  but  for  air  the  vari« 

ation  of  the  nondimensional  group  Pr.  1  c_|iA  is  minor  compared  with  variations  in 

A  P 

such  individual  properties  as  the  dynamic  viscosity  p  and  local  density  p.  This 
fact,  coupled  with  the  observation  that  only  the  temperature  insensitive  product 
Pp  centers  directly  into  the  transfoimed  compressible  boundary  layer  equatipns,  led 

Chapman  and  Rubesin  (J^)  to  the  conclusion  that  for  any  T^(x),  the  temperature  field 
within  'tee  boundary  layer  would  bear  the  same  relation  to  x  and  the  stream  function 
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T  as  in  the  low  Mach  number  casei  Consequently,  except  for  a  small  GOrrectiOO  to 
the  heat  transfer  coefficient  accounting  for  variable  p|j.,  the  procedure  outlined 
above  for  small  Mach  numbers  becomes  approximately  valid  for  arbitrary  Mach  number, 
turning  to  the  radiation  cooling  problem  at  hand  this  implies  that  all  previous 
results  for  u  =  2  (laminar  boundary  layer  flow)  can  be  extended  to  the  compressible 
case  by  (i)  replacing  the  stream  temperature  T-  wherever  it  appears  with  the  recovery 
temperature  and  (ii)  correcting  the  heat  transfer  coefficient  NUj^^(L)  for  the 
variability  of  This  procedure  should  lead  to  accurate  results  provided  the  latter 
correction  is  not  large/  However,  if  it  is  large  then  the  basic  relation  between 
the  local  heat  transfer  coefficient  and  T^(x)  must  first  be  generalized  to  account 
for  streamwise  variations  in  the  pp,  correction  (or  reference  temperature;  cf,  foot¬ 
note  below).  To  the  writer's  knowledge  this  has  never  been  Garried  throughi  althou^^ 

thhe  problem  is  far  from  intractable.* 

While  the  theoretical  justification  becomes  weaker,  similar  comments  can  be  made 
for  the  turbulent  boundary  layer  case;  i..e.  to  extend  w  =  5  results  of  Section  2,2 
to  the  compressible  case  replace  T^  everywhere  by  the  turbiUent  recovery  temperature 
and  correct  ^^^^^(l)  for  variable  properties  using  the  reference  temperature  method 

« 

subject  to  the  limitation  that  Pr^^  «  const, 

*this  correction  can  only  be  arrived  at  by  comparison  with  exact  solutions  or  experi¬ 
ment,  For  practical  applications  the  Rubesin-fckert  reference  temperature  method  has 
been  found  to  be  convenient  in  a  wide  variety  of  circumstances  (33,34), 

*for  example,  methods  similar  to  those  recently  applied  by  Mayer  (^)  could  be  used 
for  generalizing  Sq,  (11)  to  account  for  streamwise  variations  in  reference  temper¬ 
ature.  Since  the  latter  will  be  expressible  in  terms  of  9  itself,  the  remainder  of 
the  calculation  would  proceed  as  in  Section  2,2 
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(or  an  equivalent  method).  In  this  ease  it  should  be  ahtieipated  that  the  variable 

property  GorreGtion  will  quicicly  dominate  the  isolated  effeGt  of  variable  SurfaGe 

tefnperature  [or  ©n  Nu  (cf.  Fig.  4)*  This  should  allow  the  use  of  ealeulation 

methods  which  completely  ignore  the  surfaee  temperature  "history"  effect  aGcounted 

for  in  Section  2.2.  In  unusual  circumstances«  however,  the  accuracy  of  sueh  pre** 

dictions  is  questionable  (an  example  might  be  in  the  immediate  neighborhood  of  an 

abrupt  Change  in  surface  emissivity).  Here  again,  these  areas  remain  largely 
* 

unexplored. 

Summarizing  then,  in  a  chemically  inert  high  speed  stream,  radiation  losses 

will  prevent  most  of  the  surface  of  an  iitimersed,  insulating  solid  from  reaching  the 

gasdynamic  recovery  temperature  T^«  According  to  continuum  boundary  layer  theory, 
however,  the  leading  edge  region  of  sharp*nosed  bodies  will  seek  the  full  recovery 
temperature,  while  local  breakdown  of  boundary  layer  theory  and  internal  heat  con'^ 
duetion  modify  this  conclusion  somewhat  (cf.  Section  2.2)  these  effects  only  sli^tly 
increase  the  flight  Atech  number  beyond  which  localized  material  degradation  may  be 
expected.  But  for  particular  applications  in  which  localzed  failure  can  be  pre» 
vented  (say,  by  moderate  blunting)  one  has,  in  radiation  cooling,  the  possibility 
of  a  relatively  simple  and  efficient  heat  protection  scheme.  Of  course  for  long 
duration  trajectories  this  implies  the  existence  of  (i)  effective  means  for  insula* 
ting  the  radiating  vehicle  skin  fr<Hn  the  underlying  structure  and  payload;  and  (ii) 
adequate  means  for  sustaining  aerodynamic  and  aerotheimoelastic  loads  acting  directly 


an  analagous  situation  exists  in  the  case  of  chemical  surface  reactions  when  the 
catalytic  activity  undergoes  abrupt  streamwise  variations.  Several  such  eases  hove 
been  quantitstively  investigated  by  Chung,  Liu  and  Mirels  (^).  a  qualitative 
discussion  of  this  general  problem  may  be  found  in  Ref,  (27). 
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Oh  the  vehicle  skin*  Further  GOtnplicatioh  arises  from  the  fact  that  curreht  temper* 

e 

ature  limitatiohs  on  refractory  metal  coatings  renders  the  use  of  ceramic  leading 
edges  and  nose  caps  more  likely.*  As  will  be  seefli  some  of  these  ceramic  nose  and 
leading  edge  materials  can  operate  at  temperature  levels  Sufficient  to  dissociate 
oxygen  molecules  at  realistic  stagnation  pressures.  With  this  in  mind,  it  is  then 
of  interest  to  examine  the  radiation  cooling  of  solids  in  the  presence  of  thermo- 
chemical  change. 


available  coatings  to  protect  refractory  metals  from  Catastrophic  oxidation  rates 
fail  at  temperatures  well  below  the  melting  points  of  the  metals  themselves, 

*Thus,  the  case  of  preeewise  discontinuous  emissivity  alluded  to  earlier  is  of 
more  than  academic  interest. 
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3.  THERWCHEMICAL  EFFfGTS  IN  RADIATION  COOLING 


In  the  pirevidue  section  radiation  cooling  and  Gheinical  sufface  catalysis  wesre 
pfesented  as  analogous  but  separate  probleitisi  Moreoveri  the  thermal  effects  which 
aGcompany  most  Ghemioal  reactions  were  suppressed  by  considering  only  isothermal 
catalyst  surfaces  with  T  =  T.^  But  if  the  surface  reaction  is  exothermiCi  in  order 
to  prevent  increases  in  the  surface  temperature  one  must  remove  energy  continuously 
(say*  by  internal  cooling).  In  the  absenGe  of  such  cooling  the  surface  temperature 

would  rise  until,  at  each  point,  heat  loss  by  eonvection  (to  the  cooler  stream)  and 

• 

by  radiation  (e  dT  ^)  balanced  the  heat  R"Q  generated  per  unit  catalyst  area  in  the 
w  w 

chemical  surface  reaction. 


.3,1.  Generalization  of  the  Recovery _teif>perature.  Concept  (37-4Q) 


consider  first  the  extreme  case  of  negligible  radiation  loss  from  a  catalyst 
of  extremely  great  activity  (k  -►pO)  in  a  low  speed  stream.  One  then  finds  that  the 

W 

catalyst  will  settle  out  at  a  temperature  which  is,  in  general,  higher  than  T^  and 
different  from  the  mainstream  "total”  temperature 


r  -  m  Tl  ^ 


(32) 


For  lominar  boundary  layer  flow  this  catalyst  "recovery"  tempezature  is  found  to  be 


(33) 


* 

when  the  heat,  Q,  released  per  unit  mass  of  transformed  reactant  is  positive 
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where*  fey  analogy  with  the  nonreaGtiye  problem  discussed  in  Section  2.3*  the  function 
*D^^*X*^*D^  called  the  recjOjiyierv  Jaetor  . for, Jj:ee  _stream  chemiOai  energy 

Explicitly*  is  given  by 


3 

% 

(34) 


[the  last  (approximate)  equality  following  from  polhausen's  observation  (Jl)  that 
^^is^  Pr’^^  for  Pr  ^  0(1 )i]  When  the  mainstream  contains  both  kinetic  energy 
(it  u  *  per  unit  mass)  and  chemical  energy  ( o6  Q  per  unit  mass  of  mixture)*  one  then 
anticipates  the  combined  result 


* 


1 
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(35) 


which  has*  in  fact*  been  derived  by  Vaulin  (^)  for  the  special  case  of  compressible 

laminar  flow  Over  a  flat  plate.  Thus,  only  when  the  surface  temperature  T  has  the 

w 

value  T^*  given  by  Eq. (35)  does  the  net  energy  flux  to  the  plate  vanish.  If  the 
plate  is  radiating  energy  to  its  surroundings,  this  energy  must  be  derived  from  the 
main  stream.  It  follows  that  for  any  x  >  0  radiation  loss  will  cause  the  plate 


a 

The  consequences  of  these  relations  for  the  recovery  temperature  [cf.  Egs.(3i»33)] 
are  somewhat  startling  when  Pr^  >  1*  Le  >  1  (M»^)*  For  pure  gases  like  water  vapor, 
the  recovery  temperature  exceeds  the  thermodynamic  total  temperature  of  the  gas 
stream.  Similarly*  when  Le  >  1  (as  in  the  ease  of  weakly  dissociated  diatomic  gases)* 
the  steady  state  surface  temperature  exceeds  the  thermodynamic  total  temperature. 
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temperature  to  fall  below  the  generalized  recovery  temperature  T^. 

paralleling  that  given  in  Section  2i3  reveals  that  the  law  governing  this  temperature 

deeay,  T  (x)s  would  be  precisely  that  given  for  the  incompressible,  nonreactive 
w 

radiation  cooling  problem  (cf^  Section  2.2)  provided  one  replaced  by  T'_  and 
corrected  the  heat  transfer  coefficient  for  fluid  property  variations^  Thus 


=  % 


tt-j  n) 


where 


M/ 


4 


3.2  Effects  of^Finj^te  Gatalyst  Activity 


The  problem  becomes  considerably  more  complex  if  due  account  is  taken  of  the 
finite  catalytic  activity  of  real  surfacesi  Then,  even  in  the  absence  of  radiation 
loss  a  nontrivial  di stribution  of  surface  temperature  results.  in  particularj 
T  (O)  attains  the  (nonreactive)  gas  dynamic  recovery  temperature  T  and  T  (») 
attains  the  generalized  recovery  temperature  T^.  If  the  chemical  rate  constant 
has  a  temperature  dependence  of  the  two-^parameter  Arrhenius  form 


<  -  4  VKp[ ~  (38 
then,  for  a  prescribed  reaction  order,  the  surface  temperature  distribution  T^^(^) 
will  depend  upon  the  activation  energy  parameter  e/(rTj,)  as  well  as  the  catalytic 
parameter ^ ^ of  Section  2>2.  Qualitatively,  exothermic  surface  reaction  does  not 
cause  the  leading  edge  temperature  rise  above  T^  because  of  the  excellent 


* 

except  when  intemal'^heat  conduction  renders  the  surface  temperature  constant, 
thereby  allowing  isothermal  catalysis  solutions  to  be  applied  to  the  radiation 
cooling  problem  with  exothermic  surface  reaction. 

^The  rate  constant  k  appearing  in  the  definition  of  the  catalytic  parameter  ^ 
to  be  evaluated  at  the  gasdynamic  recovery  temperature  [cf.  ig.  (31)]. 


is 
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conditions  of  fteat  removal  in  this  region*  Further  downstream,  however,  the  temperature 
is  able  to  rise  above  and  this  is  accompanied  by  a  Corresponding  increase  in  the 
reaction  rate  coefficient  [cf*  Eq*  (38)1*  Very  far  downstream  the  steady  state  rate 
of  heat  evolution  (and  hence  heat  removai)  becomes  diffusion  limited  so  that  the 


surface  temperature  approaches  asymptotically*  But  this  latter  condition  would 

not  be  achieved  in  the  presence  of  radiation  heat  loss,  for  then  it  is  readily  shown 

.  # 

that  vanishes  as  X“*^in  accord  with  the  asymptotic  behaviour 


When  compared  with  the  corresponding  result  in  the  absence  of  the  chemical 
eontribution,  Eq.  (39)  reveals  that  the  addition  of  chemical  energy  to  the  free  stream 
would  increase  catalyst  surface  tenperatures  very  far  downstream  by  the  factor 


(40) 


Apart  from  mentioning  these  general  features,  we  shall  not  dwell  on  this 
coupled  radiation  cOoiing»catalysis  problem  further,^  but  instead  turn  to  a  simpler 


but  perhaps  more  practical  example  which  exhibits  similar  features.  Consider,  the 
case  of  radiation^cooled  blunt  leading  edges  or  nose  caps,  again  within  the  framewortc 


here  r  and  n  are  defined  by  Eqs.  (13)  and  (^),  respectively,  in  Eq.  (14)  is 
replaced  by  T^, 

^a  detailed  investigation  of  this  particular  problem  is,  in  fact,  nonexistent. 
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of  feOundafy  layer  theory.  Here*  in  the  usual  situation  atoffiic  oxygen  and  nitrogen 
represent  the  ehemical  energy  content  Of  the  gas  at  the  outer  edge  of  the  boundary 
layer  and  first  order  (n  =  1)  heterogeneous  reGombination  is  the  exothermiG  surfaee 
reaction  of  interest  (L2-JJi,) .  Except  for  certain  extreme  flight  conditions  or 
test  conditionsj*  atomic  oxygen  is  usually  absent  in  the  undisturbed  medium  but  is 
formed  behind  bow  shock  waves  at  a  finite  rate,  if  ^  is  the  effective  flight 
velocity  thenj  the  energy  equation  applied  between  upstream  infinity  (subscript©®) 
and  the  forward  stagnation  point  (subscript  e)  provides  the  relation 


(41) 


Thus*  for  a  given  flight  condition,  the  static  temperaturei  at  the  outer  edge 
of  the  boundary  layer  decreases  with  increasing  local  atom  concentration.  To 
understand  the  effects  Of  thermochemical  change  on  the  surface  temperatures  of  such 
radiation  cooled  solids  one  must  then  account  for  the  fact  that  the  chemical 


contribution  to  the  enthalpy  makes  its  appearance  only  at  the  expense  of  the  thermal 
contribution.  It  may  also  be  necessary  to  account  for  the  property  that  a  sufficiently 
hot  surface  can  dissociate  incident  molecules  as  well  as  recombine  incident  atoms 
(23.42) .  The  dominant  effects  of  these  phenomena  are  readily  perceived  using  the 
following  simple  model. 


3.3  Radiation  Cooled ^Stagnation  Region  wi^  Shock  Laver  Dissociation  (4) 

The  steady  state  temperature  of  a  nose  cap  will,  as  before,  be  the  temperature 
at  which  the  radiation  term  k  0~  1  balances  the  aerodynamic  heat  flux.  When 


* 

only  nonablating  surfaces  are  considered  here. 

Corresponding  to  photodissociation  or  predissociation  caused  say,  by  a  nuclear 
blast  (1^.^). 

^corresponding  to  recombination  nonequilibrium  in  shock  tunnel  nozzles  (47). 
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the  gas  at  the  outer  edge  of  the  boundary  layer  is  partially  diSSOGiated  the  input 
heat  flux  due  to  GoweGtive  proGesses  within  the  boundary  layer  may  be  represehted  as 
the  Sum  of  two  GOntri but ions,  and  Mere  represents  the  contribution  due  to 
thermochemiGal  changes  across  the  gas  phase  boundary  layer.  When  all  of  the  reaction 
ocGurs  at  the  gas/solid  interface,  each  of  these  contributions  may  be  expressed  in 
terms  of  its  nondimensional  transfer  coefficient  (Stanton  number)  and  driving  force 
as  follows  (^) 

si\- )  Q, 

where  represents  the  gas-dynamic  recovery  temperature,*  o  represents  the  mass 
fraction  of  atoms  present  in  the  partially  dissociated  mixture,  and  Q  is  the  heat  of 
atom  recombination  (assumed  constant).  The  remaining  symbols  and  subscripts  are 
defined  in  the  nomenclature. 

With  these  preliminary  statements  in  mind  the  steady  state  heat  balance  equation 
may  be  written 

Qi  (44) 

At  this  point  it  is  convenient  to  introduce  the  following  notation. 

p  ,  U5) 

7  ^  ^  0/ ^ 


(42) 

(43) 


» 

The  conditions  under  which  this  assumption  is  valid  will  be  explored  in  Section  3,4« 
The  case  of  local  thermochemical  equilibrium  yyithin  the  boundary  layer  will  be 
approximately  equivalent  to  a  special  case  ( ^  s  ©o)  treated  here, 

^At  the  forward  stagnation  point  the  kinetic  energy  term  iu  *  vanishes  so  that 

Tj  ^  Tg  regardless  of  the  magnitude  of  the  Prandti  nunber*  ®  Prj^,  Despite  the  fact 

that  p  ij  —V  P  the  product  ^\*P'U  is  finite. 

/e  -'e  c 
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(46) 


Co<^  -  o(^)  / 


(47) 


At  any  point  in  the  flow  field  the  quaifttity  T°  will  be  recioqnized  as  the  static 
temperature  which  would  be  attained  locally  if  the  kinetic  enerqy  lu^  and  the 
chemical  enetgy  a  q  were  aditbatically  converted  into  sensible  temperature  rise^ 

The  quantity  r^t  being  a  ratio  of  transfer  coefficients  (St)  for  mass  and  heatj  may 
he  regarded  as  the  recovery  factor  for  chemical  energy  and  will  be  different  from 
unity  if  the  Lewis  nuitiber  Le  for  atom  diffusion  through  the  mixture  is  different 
from  unity.  Lastlyi  the  coefficient  ^  is  a  measure  of  the  depletion  Of  atoms  across 
the  boundary  layer  due  to  chemical  change,  in  terms  of  these  quantities  Eq#  (44) 


becomes 


T 


Now  imagine  the  same  surface  element  in  a  gaseous  environment  which  does  not 
dissociate.  This  condition  will  be  designated  hereafter  by  the  subscript  Zi  for 
"inert".  In  the  "inert"  case,  if  the  flight  speed  and  other  properties  of  the  mediui 
are  unchanged,  the  surface 
equation: 


T^  j  will  satisfy  the  simplified  heat  balance 


'nf,r  =  9  (49) 

where  the  quantities  St*.  ,  C  and  T  °  will  be  approximately  equal  to  those 

Ayg  V  p  IB 

appearing  in  iq.  (46).  Identification  of  the  "total"  temperature  is,  in  fact,  a 
consequence  of  energy  conservationi  ij.e..  the  stagnation  enthalpy  of  the  gas  with 
respect  to  the  vehicle  should  be  the  same  in  both  cases.  Hereafter,  the  common 
"total"  temperature  will  be  designated  T®.  it  will  also  prove  convenient  to 
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define  the  following  dimensionless  temperature  ratios: 

ff  mr^/T^ 


i  +  -iji 


0(eQ 


T' 


€  m  T 
z 


T 


(50) 

(51) 

(52) 


and  a  GonveGtiOn-radiatiOn  pafameter  of  the  form  defined  by  Egt  (l4)f  viz* 


(T* 


(53) 


From  the  defining  relation  for  the  surface  tempefature  we  then  have  the 

familiar  n^^“degree  algebraiG  equation  for  ©j  Of) 

5te'’  »  /-  t? 

Similarly  Eq.  (48)  takes  the  more  general  form 

^  m 

whioh  reduces  to  iq*  (54)  only  when  >  1. 


When  the  convection-radiation  parameter,^  vanishes  Eq.  (54)  shows  that 
achieves  its  maximum  value  of  unity.  Thus  the  total  teiTiperature  T°  is  seen  to  be  the 
maximum  value  of  the  surface  temperature  for  motion  through  an  inert  atmosphere* 
Similarly,  when  ^’=>0  then  B  ^  J attains  the  value  ©*  in  a  dissociating 
atmosphere.  It  is  seen  from  Eqs.  (54)  and  (55)  that  the  new  temperature  variable 
T  '  i/i*  will  satisfy  the  same  algebraic  equation  as  ©j  if  the  convection-radiation 
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parameter^  is  replaced  by  ;  i.e.; 

(er'X)  -r'^  > 


1-  T 


For  PI  -  4  Fig*  8  shows  a  plot  of  the  reduced  tefflperature  %-  ©/©^t  against  the 
modified  radiation  parameter  ©J  X  •  The  cube  root  Of  this  parameter  (i.e. . 
may  be  given  an  alternative  interpretation  as  follows*  ^  itself  can  be  conisidered 
to  be  the  ratio  of  the  total  temperature  T°  to  a  characteristic  convectidn-radiation 


temperature  T^  defined  by 


(58) 


'A  4^*4  /(^^<^)  ] 

Thusi  the  parameter  i®  equivalent  to  the  grouping: 

( T'V'r^  )! i*  [rj^cj,  ~i][ 

Equations  (54)  and  (55)  reveal  that  ©  can  actually  exceed  ©j  provided  the  character*' 
istic  temperature  ratio  ©^  defined  by  Iq,  (51 )»  exceeds  unity*  This  might  be  called 
a  temperature  "overshoot"  since  the  surface  temperature  T^  would  then  actually 
exceed  that  attained  in  an  inert  atmosphere  at  the  same  velocity  (  )  and  density 

^  p  ).  Since  an  overshoot  is  impossible  when  ©^^  ^  1  the  quantity  ©^f  plays  the  role 
of  a  discriminant.  If  this  discriminant  were  rigorously  constant  then  the  extent  of 
the  temperature  overshoot  could  not  exceed  that  achieved  in  the  absence  of  radiation 
loss  (s:?->  0).  In  practice,  however,  ©^  will  itself  be  a  function  of  the  surface 
temperature  ©  and  hence>  implicitly*  a  function  of  the  convection-radiation  parameter 
^  .  Qualitatively  speaking,  the  nature  of  this  dependence  usually  makes  an  overshoot 
impossible  when  =  0  but  possible  for  large  values  of  provided  the  surface 

is  an  active  catalyst  for  atom  recombination  and  the  average  Lewis  nuaber  for  atom 
diffusion  exceeds  unity* 

This  can  be  seen  as  follows  (4).  in  the  absence  of  gas  phase  chemical  reaction, 
the  rate  of  convective  diffusion  of  atoms  to  the  gas/solid  interface  must  be 
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identically  equal  to  tlie  net  interfacial  rate  of  atom  recombination.  Owing  to  the 
possibility  Of  dissociation  at  the  radiating  surface  the  ne^  rate  of  atom  re^ 
GOmbinatiOn  must  be  written  (^*42) 

^ ■  I (t^  i  P)J  (59) 

where  the  rate  constant  k  is  related  to  the  recombioation^OOefficient  y  and 

w 

atomic  mass  m,  by  (^j42) 

)  J  •  y(% )  m 

If  the  local  atom  number  density  n^  is  replaced  by  the  atom  mass- fraction.  CC*  the 


atom  conservation  equation  at  the  gas/solid  interfaGe  becomes 

y 

^  11  '/  ^  ^  jO  J  a/  /'a  \  I  ^ 


W 


(61) 


Defining  a  non-dimensional  catalytic  parameter  solving 

for  ^  we  (A)* 


tf^) 


1  ^ 


OCeu(^) 


(62) 


a  T-  ^  ^ ^ 

where  both  71  and  l\L  (T  ;  p)  can  be  strong  functions  of  the  surface  temperaturet  T 
^  ^eq  w 


As  the  dimensionless  surface  temperature  ©  decreases  (with  increased  radiation 

loss)  it  is  seen  that  o  (©)  will  eventually  become  negligible  compared  to  o  and, 

“  '  '  -  eq  ' 


*If  the  recombination  reaction  is  truly  first  order  (n  »  l)  then  the  recombination 
coefficientF  for  a  particular  surface  will  depend  on  surface  temperature  alone. 

^effects  of  this  temperature  dependence  will  be  discussed  further  in  Section  3*5. 
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ipso  facto,  negligible  Gompared  to  unityi  In  this  extreme*  we  have  the  familiar 
result  (38.43t4Al  for  the  extent  of  recombination; 


'  1+ ^ 

On  the  other  hands  if  we  consider  the  stagnation  region  with  the  atom  conGentratiOn 


at  equilibrium  with  the  local  tempefature  it  is  apparent  that  ^  will  vanish 
when  since*  in  this  case,  Ogq(T  )  ®  Og  regardless  of  the  magnitude  of 

Thus,  we  have  a  realistic  situation  for  which,  in  the  absence  of  radiation  loss*  there 


can  be  no  temperature  over-shoot  even  for  a  "perfect"  catalyst  (i,.e, ,  even  for^»  l)» 
This  high  temperature  behaviour  of  the  extent  of  recombination*  (j)  ,  is  illustrated  in 


Figi  9  for  the  case  of  partially  dissociated 


^  at  a  pressure  of  0^1  atmi  and 


real  gas  stagnation  temperature  of  3000  K»  The  combined  effects  of  surface  reaction 
and  radiation  loss  on  the  steady  state  temperature  of  a  constant  emissivity  solid 
are  illustrated  in  fig*  10*  again  for  the  case  of  equilibrium  hydrogen  (T,  “  3000®Ki 
P  -  0»1  atm).  Here  contours  of  constant  surface  temperature  are  shown  on  the  log 


logjy  plane  together  with  a  (shaded)  domain  within  which  the  surface  temperature  is 
actually  higher  than  it  would  have  been  in  the  absence  of  dissociation.  This  overshoot 
occurs  only  for  active  catalysts  if  the  recovery  factor  r^  exceeds  unityi  and 
radiation  loss  )  is  sufficient  to  depress  the  surface  temperature  T^  to  the 
extent  that  the  corresponding  equilibrium  degree  of  dissociation  oL  {j  )  becomes 

WVgq  yij 

small  compared  to  the  atom  concentration  established  at  the  outer  edge  of  the 


While  0*is  taken  to  be  the  parameter,  it  should  be  remembered  that,  for  any 
particular  surface,  Q  will  itself  vary  as  the  surface  temperature  is  reduced. 

^Recent  transport  property  estimates  (M)  indicate  that  the  Lewis  nirnber  for  weakly 
dissociated  hydrogen  at  3000^  is  about  1.65.  This  would  imply  r^  ^  1.35, 
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boundary  layer.  On  most  of  the  log  L»  -  log^  plane  surface  temperatures  are 
considerably  lower  than  they  would  have  been  in  the  absenGe  Of  diSsOciationj  as 
refleGted  by  the  positive  slope  of  s  const.  cofitGurs  outside  the  shaded  region. 

At  point  Pi  for  examplei  the  difference  amounts  to  about  250\.  The  straight  line  of 
unit  slope  passing  through  point  P  represents  the  locus  of  operating  points  under  a 
contemplated  change  in  scale  alone,  under  such  a  Ghangei  the  positive  slope  of  t  = 
consti  contours  indicates  that  the  equilibrium  surface  temperature  will  change  by  a 
smaller  amount  than  in  the  nondissociated  case;  a  characteristic  which  will  be 
explored  more  fully  in  lection  3.4‘ 


Simple  considerations  such  as  those  presented  above  allow  one  to  single  out  the 
important  physicochemical  parameters  governing  this  class  of  radiation-cooling 
phenomena#  It  is  then  possible  to  estimate  the  magnitude  of  these  nondimensionai 
parameters  )  as  encountered  in  high  speed  motion  through  the  iarth’s 

atmosphere  by  collectively  drawing  in  the  work  of  Refs*  A  convenient 

representation  is  obtained  by  constructing  contours  (level  lines)  of  cpnstant^^^j^ 
on  the  altitude*velOCity  piane,  on  which  typical  re-entry  trajectories  Can  be  super¬ 
imposed.  Three  such  contour  maps  are  shown  in  Figs.  ll>  12>  13>  each  of  which  contains 
the  re-entry  trajectory^  of  a  typical  satellite  (Ji)  (marked  S)  and  boGst  glide 


vehicle  (^)  (marked  SGV),  The  point  on  each  trajectory  at  which  the  product 
^  d>,l%  passes  through  a  maximiBn  is  marked  with  an  open  circle 

(this  would  be  the  point  of  peak  convective  heat  input,  in  the  case  of  complete 


increasing  scale  (nose  radius)  increases  both  transfer  coefficients  St  and  Stj^  at 
equal  rates, 

^also  included  in  Fig,  11  is  a  typical  intercontinental  ballistic  missile  trajestory 
(ma^Ked  ICBM)  for  which  peak  heating  rates  are  seen  to  be  from  1  to  2  orders  of 
magnitude  greater  than  in  the  ease  of  glide  vehicle  or  satellite  re-entry. 
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atom  recombi fiat ion  with  Also  inclucied  are  dashed  GontourS  of  GOfiStant 

meximuit)  lamifiar  coftveGtive  heat  flux  (in  BTlj/ft^-seGi )  to  a  strOfigly  Gooled  surface 
of  1  foot  hose  radiusj  as  GalGulated  from  the  Detra-Hidalgo  Gorrelation  (iB). 


Figure  11  shows  the  GonvectiOfi-radiatiOfl  parameter  ^  [cf^  iq»  (53)]  for  a 

_  -  - 

nose  radius  of  1  fti  with  a  surface  emissivity  of  0*5  .  It  is  seen  that  during  re** 
entry  deGreases  from  very  large  values  [of  order  10^  at  satellite  velocity]  to 


values  of  order  10"’  at  slightly  supersofiic  velocities,  with  peak  aerodyftamie  heat 
input  OGGurrifig  at  Q(lO'i).  Figure  12,  shows  contours  of  Gonstaftt  Gstalytio 
parameter^  derived  from  the  calculations  of  Goulard  (^)  for  chemioally  frozen 

boundary  layer  flow  at  200  and  250  Kft*  (with  a  wall  temperature  of  700°k)* 
NumeriGal  values  on  each  set  of  contours  also  Gorrespond  to  a  body  with  one-foot 
nose  radius  and  the  recombination  Goeffieient  /  has  been  set  at  the  intermediate 


value*  10"^  (thus,  in  practice,  a  factor  of  10^  in  either  direction  is  possible, 
depending  upon  whether  the  nose  material  is  a  very  good  or  very  poor  atom  re-» 
combination  catalyst.)  Contrary  to  3^  ,  the  catalytic  parameter  ^  does  not 
change  very  rapidly  along  a  re-entry  trajectory,  its  absolute  value  could  be  any¬ 


where  from  10"^  to  10^  depending  on  the  magnitude  of  the  atom  recombination 
coefficient  ^  ,  Finally,  one  requires  an  estimate  of  or  at  least  its  minimuiii 
value,  since  dissoeiation  can  reduce  the  steady  state  surface  temperature  appreciably 
only  when  6^^^  is  less  than  unity  and  is  small  (4)^  Figure  13  gives  contours 


Since  3L  values  corresponding  to  any  other  nose  radius,  Rg,  or 

emissivity,  6^,  could  be  constructed  from  the  values  given  on  Fig,  11, 

*These  contours  are  marked  vy  s  const,  anticipating  the  nomenclature  to  be  adopted 
in  Section  3,4, 

*At  any  other  surface  temperature  these  contours  would  still  he  applicable,  but  to  a 

different  ,  in  accord  with  the  proportionality  k  ocyi  [cf,  Eq,  (60)}, 

w  vy 
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6f  Gonstant  1  -  q/(  ^  t®)],  where  is  the  value  of  Gomputed  behind 

normal  shock  waves  in  airi  These  values,  obtained  by  crosspiotting  Huber's  results 
for  the  temperature  ratio  aGrOss  normal  shOGks  (^),  Gan  be  regarded  as  providing  a 
lower  bound  to  the  diSGriminant  It  is  seen  that  ©«  initially  takes  on 
values  of  about  OiZ  (near  satellite  speed)  but  rises  rapidly  to  unity  near  MaGh  5» 
combining  the  magnitudes  of  and  shown  on  Figs*  li  and  13*  respeGtively, 

one  finds  that  the  actual  wall  temperature  for  the  boost-glide  vehicle  (iGV)  entry 
shown  can  be  lower  than  that  neglecting  shock  layer  dissOGiat  ion  by  as  much  as  about 
6C)0®K.  This  occurs  where  the  noncatalytic  wall  temperature  would  be  about  1400 
for  the  remaining  conditions  assuised  =  0*5,  =  l  ft)*  ®y  the  same  token  T,„ 

Can  be  larger  than  the  temperature  ealculated  neglecting  heterogeneous  atom 
recombination  by  some  60  0®K*  These  rough  estimates  suggest  that  a  low  atom 
recombination  coefficient  ^  for  high  temperature  ceramics  (or  refractory  metal 
coatings)  could  play  a  role  comparable  to  high  emissivity  in  contributing  to  the 
success  of  an  overall  design. 


3.4  influence  of -the.  Aerodvnamie  Heat  Input  Law;  -Dependence  on  Physical  gcale  (l5) 
It  has  already  been  demonstrated  that  the  steady  state  surface  temperature  of  a 
radiation  cooled  solid  should  be  less  sensitive  to  contemplated  changes  in  physical 
dimensions  than  in  the  nonreactive  case  [cf*  Section  3*3]*  This  property  can  be 
traced  back  to  the  scale  dependence  of  the  aerodynamic  heat  input  q"  since  T^  is 
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always  determined  from  an  equation  of  the  form 

-  -  *  c 

Differentiating  this  expression  with  respect  to  boundary  layer  thicknessi  o 
gives 


S' 


m 

n. 


(64) 


where  we  have  written 


9  " 


and,  as  before 


n,  m  4  ^  S  /n  6^  / c/  U  Xr 


(66) 


(2) 


The  magnitude  of  n  (i>e.  the  effect  of  the  emissivity^temperatufe  relation)  has 
already  been  discussed  [cf^  Table  l}.  The  remaining  factor  influencing 
is  therefore  m  (i.e.  the  dependence  of  input  heat  flux  q"  on  scale).  While  m  s  *1 


for  equilibrium  boundary  layer  flow  and  nonreaetive  boundary  layer  flow  alike,  in 
general  m  depends  upon  chemical  kinetic  parameters  (hence  '  In  the  cases  treated 


thus  far  gas  phase  recombination  has  been  neglected,  i.e.  the  only  chemical  kinetic 
parameter  entering  the  problem  was  (j?'  The  dependence  of  m  on  the  catalytic 


parameter  is  most  readily  demonstrated  for  the  special  case  ^gq(T^) 
then  Iq.  (63)  implies 


(67) 


« 

so  that  the  results  can  be  applied  to  laminar  or  turbulent  boundary  layer  prebiems> 
by  later  introducing  ^oc  respectively. 

P  P 
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With  some  manipulation  Eqsi  (42»53»67)  provide  the  simple  result 

(p  /l/ 


—  —  i  -i- 


where  we  have  written 


[a.^  -  de^ (r^)  IQ 

~n) 


% 


(69) 


« 


For  any  prescribed  value  of  this  chemical  enthalpy  potential  parameter  Hj  iq.  (6i) 
shows  that  m  =  -1  for  either  v  =  o  (no  recombination)  or  v  =  1  (complete  recombihatiofOl 
But  at  some  intermediate  value  of  v  the  parameter  m  attains  a  maximyiii  value}  which 
iSj  however}  less  than  aero.  Simultaineous  atom  recombination  within  the  boundary 
layer  can  cause  even  larger  departures  of  m  from  the  familiar  case  of  m  =  -i.  while 
an  adequate  set  of  such  calculations  is  not  Currently  available}  we  can}  however} 
illustrate  the  nature  of  this  effect  by  drawing  on  a  simple  model  for  which  the 
complete  nonequilibrium  solution  can  be  written  down  in  closed  form.  This  is  the 
conductivit_v_,.ee.l  1  model  (H)  depicted  in  Fig.  14*  It  is  known  that  such  models 
predict  functional  dependences  which  can  be  applied  successfully  to  many  convective 


flow  problems  in  chemical  engineeringj  indeed  this  is  the  basis  of  sO*Galied  film 
theory.  Of  course}  to  carry  this  analogy  through  to  the  computational  state  a  film 
theory^boundary  layer  "dictionary”  is  required}  i.e.  a  set  of  prescriptions  for 
translating  film  theory  parameters  into  their  corresponding  observable  quantities 
for  convective  heat  transfer.  For  stagnation  flow  this  correspondence  may  in  fact 
be  established  by  matching  film  theory  predictions  with  available  boundary  layer 


» 

this  parameter  may  be  regarded  as  the  product  of  the  recovery  factor  and  a 
"chemical  ickert  Number",  It  plays  the  same  role  in  the  theory  of  heat  transfer 
with  chemical  reaction  as  the  parameter  l^•(|  yg^)/C  ^‘(Tg  *  plays  in  non* 
reactive  compressible  boundary  layer  theory. 
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solutions  in  several  important  limiting  cases. 

TKe  Gontinuuiti  ifiodel  to  be  used  here  (cf.  Fig^  14)  is  essentially  ifoadwell's 
linearized  Gonductivity  cell  model  (53)  but  extended  to  inGlude  first  order  atom  re- 
Gombination  at  one  surfaee  (tbe  "Gold  boundary")  of  arbitrary  catalytiG  activityi 
Thus,  we  imagine  two  parallel  plates  separated  by  a  distanGe  J  ,  enGiosing  a  layer 
of  partially  dissOGiated  diatomiG  gasi  The  lOGal  atom  Goneentration  ^  is 
postulated  to  attain  its  eguilibrium  value  at  the  hot  plate  temperature  andi  within 
the  cellj  thermal  diffusion  and  other  seoondary  transport  meGhanisms  are  assuiiied 
negligiblei  For  this  problem  three  important  nondimensional  parameters  emerge, 
written  G,  Wj  H.  The  parameters  6  and  M  are  preeisely  those  appearing  in  Sroadwell's 
study,  the  first  Of  these  being  a  gas„Dhase  reGombination,  parameter  interpretabie  as 
the  ratio  between  the  charaGteristiG  diffusion  time  ^ across  the  film  to  the 
gas  phase  recombination  relaxation  time  .  The  parameter  |i,  already  defined  by 
Iq.  (69)j  is  seeh  to  be  a  measure  of  the  maximuin  possible  ehemical  enthalpy  change 
across  the  cell,  as  compared  to  the  sensible  (frozen)  enthalpy  change 

CnciTl)  vQ* 

<ih,  correspondihg  to  the  imposed  temperature  difference  AT  i  T  -  T  ,  The 

4  ■  6  W 

additional  kinetic  parameter,  ,  defined  by 

flT  m  S'/" T>,^  (70) 

t 

is  introdueed  by  the  finite  rate  atom  recombination  at  the  colder  surface,  '  The 
special  ease  71^=  0  eorresponds  to  a  completely  noncatalytic  surface  (as  treated  by 


A  nonequilibrium  Couette  flow  model  of  the  type  investigated  in  Refs.  53-55_  would 
be  required  for  loeations  on  the  body  at  which  viscous  dissipation  must  be  accounted 
for. 

^making  use  of  the  film  theory-boundary  layer  dietiGnary,  this  parameter  will  be 
recognized  as  identical  to  ^  ,  This  identity  should  be  kept  in  mind  during  the 
following  discussion. 
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Sroadwell)  whereas  ^  Gorresponds  to  the  diffusion  Gontfolled  ("perfeGtly" 
catalyticj  case.  The  behavior  of  the  Nusselt  number  Isiu  (based  on  the 

temperature  differenGe  ^T  aGrOss  the  film  ©f  thiG'kness  )  is  shown  in  Fig.  15  for 
the  speGial  case  H  =  9,  that  ^  approaches  1  +  H  for  either 

large  G  w  lar&e  since  the  extent  of  recombination  e  is  complete  in  either  case. 
On  the  Other  handi  for  very  atall  values  of  Q  and  t£,  m  for  o  the  heat  transfep 
coeffiGient  approaches  its  nonreactive  value*  unity. 

Contours  of  Gonstant  ^  have  already  been  displayed  on  the  altitude-velocity 
plane  (cf.  Fig.  11 ).  in  a  similar  way*  by  drawing  on  the  work  of  Refs.  52.57 
contours  of  constant  G  and  H  may  be  constructed  to  provide  an  overall  picture  of 
the  nonequilibrium  flight  regimes  to  be  expected  in  the  earth's  atmosphere.  Figs. 

16  and  17  display  this  data  quantitatively  in  the  range:  10  <  30Kft./see»; 

100  <  h  <  300  Kft,  Fig.  16  is  obtained  from  a  cross  plot  of  Whalen's  calculations 
(f2)  of  the  gas  phase  recombination  parameter*  and  applies  tu  a  body  of  1  ft.  nose 
radius.  Numerical  values  indicate  that  lifting  vehicles  at  very  high  altitudes  will 
most  likely  experience  Chemically  froreh  flow  in  the  boundary  layer  (G  ^  0)*  or  the 
incipient  effects  of  gas  phase  recomfaihationj  Istimates  of  the  chemical  enthalpy 
potential  parameter  H  shown  in  Fig.  l7  were  derived  from  cross  plots  Of  Grier  and 
Sand's  calculations  (^)  of  the  maximum  possible  heat  flux  reductipn  associated  with 
a  completely  noncatalytic  (  V  ^  o)  wall.^  At  the  point  of  peak  heat  transfer  to  the 


outside  of  this  range  contours  have  been  sketched  in  to  show  only  qualitative 
behavior. 

^a  tendency  which  is  accentuated  at  high  surface  temperatures*  (M)>  such  as  are 
likely  to  prevail  in  radiation  cooling  applications. 

^numerical  values  in  this  case  pertain  to  the  choice  T  ^  1000%  but  H  could  be 
larger  in  practice  if  high  surface  temperatures  are  considered;  e.o.  order  10  for 
high  temperature  ceramics. 
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boost-glide  vehicle  (EGV)  the  heat  flux  oouid  be  reduced  by  almost  a  factor  Of  4. 
With  these  magnitudes  in  mind  the  effect  Of  the  Ghemicai  kinetic  parameters  G  and 
W  on  the  scale  dependence  parameter  m  can  be  estimated  using  the  conductivity  cell 
model,  for  any  prescribed  value  Of  Hi 


consider  a  change  in  scale,  or  in  the  case  of  the  conductivity  cell  model,  a 
change  in  the  film  thickness  Then,  from  the  definition  of  the  heat  transfer 
coefficient  Nu, 


V  2>hi  {  ^/n 

lir 


m) 


where  the  geometric  interpretation  of  the  partial  derivatives  is  clear  with 
reference  to  Fig.  15.  Evaluation  of  these  partial  derivatives  leads  to  results  of 
the  type  sketched  in  Fig.  18,  which  shows  contours  of  constant  in  projected  Oh  the 
log  W  -  log  G  plane^  for  H  =  9.  Similar  to  the  case  of  the  log  \o%^  plane 
discussed  in  Section  3«3,  it  is  instructive  to  consider  the  dependence  of  g  along 
a  straight  line  of  constant  W//g  passing  through  some  salient  operating  point  (g, 
Under  a  contemplated  change  in  the  film  thickness*  (  ur  )  both  kinetic  parameters 


assuming  ail  other  quantities  (in  particular  H)  remain  substantially  unchanged. 

This  ceases  to  be  true  in  the  flight  regimes  where  dissociative  relaxation  in  the 
shock  layer  becomes  important. 

^these  contours  are  also  included  on  the  isometric  drawing  of  the  Nu  (G,  Wj  9) 
surface  shown  in  Fig.  15« 

*for  a  conductivity  cell  iT  may  be  identified  with  the  volume/active  surface  ratio. 
Thus,  in  considering  the  effects  of  a  change  in^  ,  one  is,  in  a  sense  inquiring  into 
the  behavior  of  a  "reactor”  under  a  change  in  surfaee«-to-volume  ratio. 
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G  and  W  would  change  but  in  sucib  a  way  that  the  new  operating  point  (G'W )  would  lie 
along  this  line.  If  ^  were  smaller  than  I'  then  q"  at  the  new  operating  point 
(^'W' )  would  indeed  be  larger  than  at  (G,w),  however  the  heat  flux  sensitivity  to 
the  oontemplated  ohange  in  ^  will  be  considerably  less  than  proportional  to  ^ 
over  an  intermediate  portion  of  this  line  (increasing  toward  ^  only  for  sufficiently 
small  i  ).  Using  iq.  (65))  similar  statements  are  seen  to  be  valid  for  the  scale 
dependence  of  the  steady  state  surface  temperature  for  radiation  cooled  leading 
edge  elements. 

Before  concluding  with  a  discussion  of  the  effects  of  surface  temperature  on 
the  rate  of  aerodynamic  heat  inputs  several  comitents  should  be  made  concerning  implicit 
limitations  on  the  application  of  the  altitude- velocity  maps  presented  in  Figs,  lis 
ISs  13)  16,  17.  First)  altitude  limitations  are  set  by  the  breakdown  of  boundary 
layer  theory  (h '2^250  Kft*  for  1  ft.)  and  the  process  of  shock  layer  dis¬ 
sociative  relaxation  (which,  in  the  re-entry  vehicle  velocity  range  may  be  important 
as  iovN  as  200  Kft.J  of.  Refs.  M,  Second,  radiation  heat  transfer  t£  the 
surface  will  begin  to  contribute  significantly  (  >  10  percent)  to  the  total  heat 
flux  in  the  domain  roughly  charaeterized  by  q"^  lO^  BTU/ft.*  »  see.  (60).  Since 
solid  surfaces  would  have  to  achieve  temperatures  of  at  least  3750°K  to  re-radiate 
a  heat  flux  of  this  magnitude,  this  justifies  neglect  of  gas  cap  radiation  for 
most  applications  involving  cooling  by  radiation  alone.  Finally,  at  super- 
satellite  velocities  (where  cooling  solely  by  radiation  is  not  feasible  except  at 
extreme  altitudes)  ionization  phenomena  must  be  accounted  for  in  calculating  both 
the  conveetive  and  radiative  heat  flux  (6l.  62). 

3.5  Influence  of  the  Aerodynamic  Heat  Ipput  Law;  Dependence  on  Surface  Temperature 
ingygl 

A  discussion  of  heat  transfer  with  chemical  reaction  would  not  be  complete 
without  some  mention  of  the  very  dramatic  effects  that  can  be  produced  by  changes  in 

surface  temperature*  For  example,  the  input  heat  flux  can  have  the  same  value  at 

*  - ^ 

The  scale  dependence  of  surface  temperature  is  of  interest  in  the  pyrometry  of 
flame  gases  since  extrapolation  procedures  (to  zero  diameter)  are  used  to  obtain 
the  true  flame  temperature  (6).  41 
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three  distinct  surface  temperatures,  as  pcinted  out  in  Ref.  14*  As  a  cordiiary  it 
should  then  be  possible  for  three  steady  state  temperatures  to  exist  for  a 


radiation  cooled  solid  in  a  given  environment,  this  is  attributable  to  the  fact 
that  each  of  the  contributions  to  the  total  heat  flux  q"  =  "t"  q"p  has  a  distinct 

temperature  dependence,  and  the  chemical  Gontribution  (6'’,^)  can  actually  increase 
with  increasing  surface  temperature  owing  to  a  strong  temperature  dependence  of  the 
recombination  coefficient,  lehavior  of  this  type  is  illustrated  in  Fig.  19  for  the 
case  of  stagnation  point  heat  transfer  to  an  axisymetric  catalytic  solid  when  the 
degree  of  dissociation  at  the  outer  edge  of  the  boundary  layer  is  the  equilibriuiij 
Value  corresponding  to  the  temperature  T^.  As  discussed  in  Section  3.3  this  latter 
condition  implies  that  both  qj^  and  q''^  will  vanish  when  the  surface  reaches  the 
temperature  T^,  regardless  of  its  catalytic  activity  (i.e.  for  any  ^  ).  For  all 
surface  temperatures  lower  than  T^  the  heat  flux  to  a  perfect  catalyst  l)  is 


always  larger  than  the  heat  flux  to  a  Gompletely  inactive  surface  ( =  0)  in  the 
absence  of  gas  phase  recombination.  Owing  to  the  choice  of  ordinatej  this 
difference  is  represented  by  the  vertical  distance  which  attains  a  constant 

maximum  value  for  surface  temperatures  below  the  dissociation  "threshhpld".  The 
heat  flux  in  the  absence  of  any  recombination  (  ^  s  O)  increases  almost  linearly  as 
the  surface  temperature  decreases,  as  represented  by  the  growth  of  the  vertical 
distance  ^Ih^,  Now  consider  the  heat  flux  to  a  catalyst  for  which  the  rate  eonstant 
increases  rapidly  with  surface  temperature  (ef.  dashed  curve.  Fig,  19).  At  low 


surface  temperatures  k  is  so  small  that  the  heat  flux  is  shown  near  the  nonreactive 

w 


Moore  and  Pailone  (41)  have  provided  an  exact  solution  (h  s  200  Kft.,  AJxd  ^ 

=  Q»  Rg  =  1  ft.)  displaying  precisely  the  limiting  behavior  shown  in  Fig.  19, 
For  this  case  Tg  (^v6750CK)  is  far  beyond  the  destruction  temperature  of  known 
materials.  Interestingly  enough,  for  surface  temperatures  in  the  range  2000  « 
3000%,  their  solution  exhibits  nonequilibriym  properties  consistent  with  the 
contours  presented  in  Figs.  16  and  17, 
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value  (Gorresponding  to  p^iO).  At  larger  temperatures  the  extent  of  reGombinatiGn 
rises  rapidly,  possibly  attaining  a  value  close  to  unity  for  temperatures  below 
The  result  is  a  heat  input  curve  exhibiting  two  extrema  and  a  region  of  positive 
Slope  on  the  '^"‘“T^  plane,  'Following  the  arguments  of  ftef,  1^,  if  this  were  the  Gase 
a  range  of  intermediate  surface  temperatures  Gould  be  statiGally  unstable,  depending 
on  the  nature  of  the  heat  ioss*temperature  relation.  This  unstable  region  would 
certainly  be  contained  within  the  two  temperatures  at  which  the  g"(T^)  relation 
exhibits  its  extrema,  as  has  been  pointed  out  by  Carson  (^),  Indeed,  using  the 
general  correlation  equation  of  Ref.  ^  Carson  has  displayed  a  specific  case  in 
which  such  a  region  would  be  expected  to  exist.  This  is  shown  in  Fig.  20,  applicable 
to  a  LiCl  Surface  of  nose  radius  1.25  ft.  coated  with  lithiuin  chloride,  operating  at 
h  =  250  Kft.,  15  Kft./sec.  In  this  case  it  can  be  verified  that  the  non*» 


monotonic  behavior  occurs  for  temperatures  at  which  radiation  loss  alone  eoyld  not 
possibly  cope  with  the  input  heat  flux-  But  for  a  radiation  cooled  solid  in  a  fixed 
aerodynamic  environihent  two  stable  steady  state  surface  temperatures  are  possible 
only  if  at  some  intermediate  temperature  level 


>  ^ 


It  can  be  shown  from  Iqs.  (4^?, 43,63)  that  this  is  possible  only  if 


t 


(73) 


>  4n^ 


iv 


(74) 


tt 

however,  radiation  loss  would  be  a  contributer  to  whatever  heat  protection  system 
was  adopted. 

^Interestingly  enough,  the  recombination  rate  data  used  by  Carson  indicates  that 
this  necessary  condition  would  be  met  at  T  540lc if  LiCl  were  characterised  by 
n  ^  3.5  or  less  (cf.  Table  l). 
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where  l/ft  describes  the  temperature  dependence  of  the  rate  constant  in  accord 
with  iq*  (3S)  and  n  (-  4  for  conistant  emissivity  surfaces)  is  again  given  by  14*  (2)* 
TherefOrei  if  Over  the  entire  temperature  range  below  the  melting  point  the  kinetic 
parameter  i/RT^  never  exceeds  4n>  mulitple  solutions  and  hysteresis  phenomena 
(cf*  Ref*  14)  can  be  ruled  out*  When  the  exothermic  chemicai  process  of  interest 
is  surface  catalyzed  atom  recombination>  the  presence  of  the  factor  n  (in  place  of 
unity)  makes  this  class  of  phenomena  less  likely  for  radiation  cooling  than  for 
convection  cooling,  in  the  latter  case  it  should  be  remarked  that  the  phenomenon 
is  quite  commonj  in  fact  Lavrovskaya  and  voevodskii  (65)  have  used  it  to 
experimentally  determine  i  and  e/(RT)  for  oxygen  and  hydrogen  recombination  on 
several  surfaces* 


unless  radiation  loss  contributes  to  another  loss  mechanism  which  beccmes  less 
effective  at  higher  surface  temperatures. 
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4*  OONCLUDING  REMARKS 

iln  parts  2  and  3  we  nave  emphasized  aspects  of  radiation  cooling  which  have 

received  relatively  little  attention  but  which  may  in  practice  play  a  very  im,portant 

role  in  determining  steady  state  surface  temperatureSi  While  the  treatment  has  been 

simplified  so  as  to  illustrate  the  major  effects  in  most  particular  cases  more 

* 

refined  computational  methods  can  be  brought  to  bear  to  increase  accuracyi  But, 
particularly  in  the  area  of  chemical  nonequi librium  effects,  where  accurate 
computations  may  be  very  time-consuming  (and,  hencei  expensive),  an  overall  picture 
of  the  regimes  in  which  large  effects  are  possible  is  of  interest  so  that  the  heavy 
artillery  can  at  least  be  trained  Oh  the  proper  targets. 

The  analysis  also  displays  the  important  similitude  parameters  and  input  data 
requirements.  In  the  area  of  materials  properties  it  is  clear  that  in  order  to  make 
accurate  surface  temperature  predictions  high  temperature  data  is  urgently  needed  on 
(i)  emissivity,  (ii)  thermal  conductivity,  (iii)  recombination  coefficients.  Long 
extrapolations  which  are  not  guided  by  sound  theory  can  prove  quite  inaccurate  in 
each  of  these  categories.  To  illustrate  this  in  case  (iii)  we  have  only  to  look  at 
existing  recombination  rate  data  on  quartz  surfaces  (66.67),  where  marked  deviations 
from  simple  Arrhenius  behavior  is  evident.  For  nitrogen  atoms  On  quartz  at 
temperatures  near  1200°K  the  rate  constant  k  appears  to  obey  a  relation  of  the 

form  of  Eq.  (3$)  but  with  Aqc  T  i.  This  implies  the  existence  of  a  temperature  at 
which  /  would  maximize;  a  phenomenon  potentially  important  for  high  temperature 
ceramic  materials,^ 


An  obvious  improvement  would  result,  for  instance,  from  use  of  enthalpy  instead  of 
temperature  as  the  dependent  variable  in  the  convection  part  of  the  problem. 

^For  nitrogen  atoms  on  quartz  (800°-1200%)  data  (67)  indicate  that  i  would  attain  a 
maximuiTi  value  somewhat  less  than  2.5  x  1U~3  at  an  extrapolated  temperature  of 
2300^  (which  is  well  above  its  melting  point). 
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In  the  area  of  §as  phase  properties  the  demands  are  of  course  those  already 
familiar  in  aerodynamic  heating  calculations;  i. e .  better  information  on  transport 
properties  (of  multicomponent  gases)  and  chemical  kinetic  constantsj  particularly  at 
high  temperatures^  Here  even  interpolation  is  hazardous»  as  illustrated  by  the 
intriguing  result  (^)  that  the  temperature  dependent  cross  sections  for  the  H-Ha 
interaction  are  probably  not  intermediate  between  those  of  H-H  and 

Finally,  mentiofi  should  be  made  of  problem  areas  in  transport  theory  implicitly 
or  explicitly  suggested  by  earlier  examples.  These  would  include  generalized 
predictions  of  the  effects  (on  surface  temperature)  of  (i)  mutual  radiation  transfer, 
as  encountered  in  external  flow  over  locally  concave  bodies  or  for  internal  flows 
(ii)  finite  material  thermal  conductivity  [cf.  Sections  2.2  and  Ref.  M]  and  finite 
thermal  diffusivity  (for  analyses  of  transient  behavior)  (iii)  simultaneous  mass 
transfer,  e.,a.  as  encountered  in  using  charring  ablators  (iv)  piecewise  discontinuous 
emissivity  [cf.  Section  2.3]  (v)  departures  from  continuuiti  gas  dynamic  behavior  at 
high  altitudes  (vi)  homogeneous  and  heterogeneous  chemical  kinetic  phenomena  [of. 
Sections  3.2-3* 5]  and  (vii)  aerothermoelastic  coupling  phenomena  (69). 
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the  result  that  aerodynamic  heat  flux  is  approximately  dependent  only  on  enthalpy 
difference  (across  the  boundary  layer)  has  led,  with  distressing  frequency,  to  the 
erroneous  impression  that  chemical  kineties  dp  not  exert  a  strong  influence  on  heat 
transfer.  The  fact  is  that  chemical  reaction  rates  play  a  major  role  in  establishing 
the  magnitude  of  Ah  in  most  hypersonic  applications.  Fig.  17  indicates  that  local 
heat  fluxes  can  easily  vary  by  a  factor  of  depending  upon  the  extent  of  re¬ 
combination. 
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NOt^NOLATURi 

pre-exponehtial  faGtor  in  Arrhenius  expression;  iq.  (3^) 
nondimensional  catalysis  parameterj  iq>  (27) 
specif ic  heat  of  mixture  (frozen) 
atom-moleGule  binary  diffusion  coefficient 
activation  energy;  iq.  (38) 

(  S  V0i  2)/^i  99s  phase  recoiabination  rate  parameter 

D  chemi  eq'  f 
enthalpy  of  gas  mixture 

altitude 

Boltzmann  constant 

first  order  rate  constant  for  heterogeneous  recombination 

Lewis  number  i  Dij/t  X/(  |>Cp)3 

total  length  along  surface 

exponent  defined  by  tq.  (66) 

atomic  mass 

Mach  number 

Nusselt  number 

true  reaction  order*  Eq.  (5);  or  exponent  defined  by  Eq.  (2) 
atom  number  density 


heat  of  recombination 

recovery  factor  for  free  stream  kinetic  energy;  Eq«  (31 ) 
recovery  factor  for  chemical  energy;  iq.  (33) 
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R 

B  nose  radius  of  body 

R  universal  gas  constant 

^  nondimensional  GonVectiOn-radiation  paraineter!  Ig.  (l4) 

local  rate  Of  reaGtiOn 
St  Stanton  nuiTiber 

t  wall  thickness 


T  absolute  temperature 

T®  total  temperature  defined  by  Iqs»  (29»32j45) 

Tg  Gharacteristic  convection-radiation  temperature;  Sq^  (5?) 

u  component  of  local  gas  velocity  parallel  to  surface  element 

\j^  velocity  of  vehicle  relative  to  undisturbed  atmosphere 

W  catalytic  parameter  defined  by  |q»  (70) 

X  distance  along  solid  surface 


y  distance  normal  to  solid  surface 

z  nondimensional  coordinate  defined  by  |g.  (15)>  cf.  Table  2 
oL  local  mass  fraction  of  reactant  (atoms) 


i 

J 

e 

w 

n 

n 

5- 

e 

©a 

X 


X' 


inviscid  velocity  gradient  at  nosei  of.  Fig.  19 
recombination  eoefficientj  Eq.  (60) 
plate  spacing,  or  boundary  layer  thickness 
total  hemispheric  emissivity  of  surface 
dummy  variable 

normalized  local  reaction  rate  or  radiative  transfer  rate;  ef.  Table  2 

catalyst  or  radiation  effectiveness  factor;  cf.  Table  2 

local  driving  force  for  diffusion  or  heat  transfer 

normalized  local  surface  teimperature 

discriminant  defined  by  iq.  (51 ) 

thermal  conductivity  of  mixture  (frozen) 

thermal  conductivity  of  solid 
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dynamic  viscosity  of  mixture 
kinematic  visGosity  of  mixture 
absolute  density  of  mixture 
Stefan-ioitemann  radiation  constant 
recombination  relaxation  time  (^)  in  gas  pbase 
©/e^^^  Eq*  (56) 

extent  of  reGombination;  Iq.  (47) 
local  transfer  rate  coefficienti  cfa  Table  2 
global  transfer  rate  coefficient;  cf.  Table  2 
Von  Mises  stream  function 

parameter;  s  2  for  laminar  boundary  layer  flow;  ^ 
boundary  layer  flow 

Subscriots 

chem  Ghemical  contribution 

D  pertaining  to  diffusion 

e  at  Outer  edge  of  boundary  layer 

eq  pertaining  to  local  thermochemical  equilibrium 

f  chemically  frozen 

G  at  constant  G 

I  inert 

iso  isothermal 

max  maximujTi 

min  minimum 

r  pertaining  to  recovery  condition;  Eqs*  (31 ) 

$  behind  shock 

w  at  the  wall  (gas/solid  interface) 

W  at  constant  W 


5  for  turbulent 
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pertaining  to  moleGular  conduction 

1  atoms 

2  molecules 

oo  st  upstream  infinity 

Operators 

Change  in 

In  natural  logarithm 

0(  )  order  of  magnitude  of 
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TABLE  I 

APPROXIMATE  tEMPBRATURE  DEPENDENCE  OF  TOTAL  RADIATION 
FLUX  FROM  HiATED  SOLID  SURFACES 


Surface 

Temperature  Range  (°K) 

n 

Ref, 

opaque  quartz  (vitreosil) 

530-1090 

3.5 

20 

alumina  ceramiG 

530-1090 

2.5 

20 

thoria  ce ramie 

530-1090 

2.3 

20 

magnesia  ceramic 

530-1090 

2.5 

20 

quartz  glass  (2  mm) 

530-1090 

2.9 

20 

iron  (micrometeor) 

1500-3000 

4.65 

21 

iron 

700-1300 

5.55 

22 

ni  chrome 

325-1310 

4.1 

22 

polished  tungsten 

420-530 

5.4 

20 

polished  magnesium 

420-530 

5.1 

20 

mo  lybdenum 

420^530 

5.2 

20 

96?  pure  polished 

aluminum^ 

420-530 

5.4 

20 

nickel 

463-1280 

4.65 

22 

pure  polished  silver 

420-530 

4.0 

20 

silver 

610-960 

4.1 

22 

pure  polished  platinum 

420-^30 

4.9 

20 

platinwi 

640-1150 

5.0 

22 

Ref.  20  also  cites  data 

for  aluminum  which  suggests 

n  «  4.9  in  this 

same  temperature  range. 
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TABLi^a 

ANALOGY  BtTWilN  SURF AGE  CATALYSIS  AND  RADI ATI ON  COOLING  FOR 
THE  FLOW  OF  A  SLIGHTLY  VISCOUS  FLUID  OVER  A  FLAT  PLATE 


Quantity  Surface  Catalysis  Radiation  Cooling 


Unknown 


OC  = 


e  = 


Independent  Variable  z  =  C  *  (x/l 


Parameters 


Local  Coefficients 


n  s  reaction  order 


*  Actual  Local  Rate 
k  (pa^)'’ 


z  =lR-(x/L) 


n  s  4  +  d(ln  e  )/d(ln  T) 

W 


3R,« 


-  Actua 1 _Loca 1 _  Rate 


w'  e 


0  = 


Actual  Local  Rate 

Convective  Rate  (Y  -*  0 
w 


0  5 


Actual  Local  Rate 
Conv  ective  Ra te ,  T^^  o 


integrated  Coefficients  r]  = 

k 


^  Actual  Overfall  Rate 

e  (T  )*aT  *L 
w'  e'  e  ^ 


.  ^  1 1  n  X  7.  -  Actual  Overfall  Rate 

Integrated  Coefficients  »  =  c^^ecti^^Rate!o^^-^^  ^  ddWectiW  RaterT^ o 


20001 


(after  Naysmith  and'  Woodliey)) 


Fig.  4  Effect  of  a  noideothennal  surface  on  the  magnitude  of  the 
local  heat  transfer  coefficient;  boundary  layer  flow 


NORMALIZED  LOCAL  RADIATION, 


Fig*  6  Normi^ized  distribution  of  mdlative  flux 
for  constant  enaitsivity  flat  plates 


for  flat  plates  with  powers  law 
emlsaivity-temperature  relation 


»n  ol 


<c  aQO 


I  10  100 

RADIATION  PARAMETER ,  ^  •» 


Fig*  10  I^eii  of  conststnt  surface  temperature 

for  raiiiation  cooled  catalysts  of 

arbitrary  activity  in  hydrogen 

(p  =  0, 1  atm,  T  =  3000^) 
c 


AUITUDC,  fcilofeet 


ALTHiTUiDE ,  kiilo  f eel 


Fig,  13  Qualitative  behavior  of  6*  ,  on  the 

■ ,  min 

altitude* velocity  plane 
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CATALYTIC  COLD  PLATi  (W) 


CONDUCTIVITY  CELL 


Fig.  14  Geometry  nomenclature; 
thermal  conductivity  cell 
(after  Broadwell,  1958) 
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Fig.  18  I^cii  of  epnstauit  m  on  the  W-G  plane  (H  "  9) 
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Fig.  20  Deipendence  of  stagnation  point  heat  Oux  to  a  Li  Cl  snrfaee 
in  the  iregiane  of  chemically  frozen,  axisymmetrie  laminar 

hovindary  layer  flow 


